Engineering entanglement for metrology with rotating matter waves 



cn 

o 

(N 

Oh 
< 

(N 



Oh 
I 



> 

00 

m 

o 



X 

c3 



L.M. Rico-Guticrrcz, T.P. Spiller, and J. A. Dunningham 
School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom 

Entangled states of rotating, trapped ultracold bosons form a very promising scenario for quantum 
metrology. In order to employ such states for metrology, it is vital to understand their detailed 
form and the enhanced accuracy with which they could measure phase, in this case generated 
through rotation. In this work we study the rotation of ultracold bosons in an asymmetric trapping 
potential beyond the lowest Landau level (LLL) approximation. We demonstrate that whilst the 
LLL can identify reasonably the critical frequency for a quantum phase transition and entangled 
state generation, it is vital to go beyond the LLL to identify the details of the state and quantify 
the quantum Fisher information (which bounds the accuracy of the phase measurement). We thus 
identify a new parameter regime for useful entangled state generation, amenable to experimental 
investigation. 



I. INTRODUCTION 

Quantum metrology is rapidly emerging as an exciting 
and feasible new technology based on quantum entangle- 
ment [1-8]. Great strides have been made in overcom- 
ing a range of practical issues that have hampered its 
implementation and proof of principle experiments have 
already been performed both with light [9, 10] and atoms 
[11] although real- world applications are yet to be real- 
ized [12]. One class of systems that has been attracting 
particular attention is atomic Bose-Einstein condensates 
(BECs) trapped in rotating potentials [13]. These are ap- 
pealing because they are within reach of current experi- 
ments [15-17], provide a conceptually simple way of gen- 
erating many-body entanglement, and have the prospect 
of leading to the development of ultra-precise gyroscopes 
[18]. 

A particularly promising system for achieving these 
goals is a BEG in a quasi two-dimensional trap with a 
weak anisotropic stirring potential. As the rotation fre- 
quency is slowly ramped up from zero, the ground state 
passes a critical frequency 0,^ signalling the emergence of 
the first vortex. At this point, the system undergoes tur- 
bulent symmetry breaking that heralds a quantum phase 
transition and the many-body state is a strongly corre- 
lated entangled state involving two macroscopically oc- 
cupied modes [19, 24]. Such systems may be amenable 
to experimentally generating entangled states for use in 
quantum metrology schemes. 

In this paper, we present a detailed calculation to as- 
sess the usefulness of the entangled states that can be cre- 
ated. We begin with a LLL calculation. This allows for 
a simplified and more computationally tractable descrip- 
tion of the system but requires that the particle density 
is low, the interaction strength is small, or the rotation 
frequency is very close to the harmonic trap frequency. 
We find that this approximation accurately predicts the 
critical frequency ilc to within a few percent. This could 
lead us to believe that it provides a good description of 
the system. However a deeper investigation that goes 
beyond LLL reveals that this is not true in general. We 
find, for example, that the 'true' many-body state at J7c 



and the quantum Fisher information, Fq, which quanti- 
fies how well a state performs in a quantum metrology 
scheme, are significantly different from the LLL predic- 
tions. 

So, while the LLL approximation may be sufficient for 
some purposes such as studying symmetry breaking [19], 
a more detailed calculation that includes a larger basis is 
vital if we are interested in the precise details of the en- 
tangled state or if we want to use it in applications. This 
issues a warning about the use of the LLL approxima- 
tion, particularly as a guide for designing experiments. 
A study of the anisotropic trap with our more detailed 
calculation reveals a rich system with interesting oppor- 
tunities for engineering different entangled states. In par- 
ticular, we are able to identify a new parameter regime 
for generating entanglement with a form very similar to 
a bat state [1, 20]. Bat states are known to be well- 
suited to metrology because they combine high measure- 
ment precision with robustness to loss. This gives them 
a distinct advantage over NOON states, which are the 
archetypal states for metrology. NOON states have the 
form (|iV)|0) -(- |0)|7V))/%/2, where the two kets in each 
term represent the number of particles in two different 
states, and are very fragile because as soon as there is 
any information about which state any one particle is in, 
the superposition is destroyed. 



II. MODEL 

The physical system we consider consists of a meso- 
scopic sample of N bosonic atoms of mass M in an ax- 
ially symmetric harmonic potential, with frequency a;j_ 
in the xy plane and u^ in the z axis, interacting through 
hard-core-type elastic collisions. We take hu:^ to be very 
large compared to the interaction energy so that the dy- 
namics along the z axis are frozen, i.e. all particles oc- 
cupy the lowest axial energy level, thus rendering the gas 
effectively two-dimensional at sufficiently low tempera- 
tures. The trapped gas is rotated at angular frequency 
Q. around the z axis with the aid of an external poten- 
tial which in the rotating frame appears as an anisotropic 



quadratic potential in the xy plane. The Hamiltonian in the rotating frame is given by 
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where A{<^ 1) quantifies the degree of anisotropy, which 
is treated perturbatively in the calculations. Here, Lzi is 
the angular momentum component in the z direction of 
the «-th atom and the ^flLz term transforms the Hamil- 
tonian to the rotating frame, where L^ is the total an- 
gular momentum of the condensate. Finally, g is the di- 
mensionless interaction coupling constant which quanti- 
fies the strength of two-particle interactions and is related 
to the 3D scattering length a hy g = vStto/A^, where 
A^ == y/h/Muj,. 

In the absence of anisotropy {A = 0), the Hamilto- 
nian is exactly diagonalizable in blocks of definite total 
angular momentum by choosing a Fock basis of the form 
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where dk creates a boson in the single-particle state 
k and A^k is the occupation number of level k. 
Here, the index k represents a pair of quantum 
numbers (n, m), the principal quantum number and 
the projection of the angular momentum respectively 
[21]. For this basis, in the strictly noninteracting 
case, the Hamiltonian is diagonal with eigenvalues 

E/Tiu!± — 2 J2i "-« + X^i l"^il ~ ^^ + ^j where L is the 
total angular momentum of the system in units of h and 
n is the rotation frequency in units oiuj±. Therefore, in 
the case of independent bosons, when Q sa ujj_ the spec- 
trum consists of highly degenerate levels called Landau 
levels that are separated by an energy gap of 2huj± . 

The LLL approximation corresponds to considering 
only the states with rn = and rrii > 0. This greatly 
reduces the Hilbert space dimension, allowing the simula- 
tion of larger numbers of atoms, as well as being very con- 
venient for computational implementations. The LLL ap- 
proximation is valid when the chemical potential is much 
smaller compared to the level spacing, i.e. /i <^ 2huj±. In 
the limit of fast rotations, fl — >■ u)±, the chemical poten- 
tial is just the mean-field interaction energy per particle, 
i.e. gn, where n is the mean density of the condensate 
and g is the interaction strength related to the dimension- 
less coupling constant hy g = h g/M. Therefore, in this 
limit, the standard validity criterion is that gn <^ 2hLUj_. 

In the present study, we will investigate the effect of 
including more Landau levels in the calculations, when 
doing so, we define the nLL-th Landau level approxima- 
tion as that one obtained by including only basis states 

where J2i=i ['^i + (l"^j| ~ "^i)/2] = J^ll — 1, as consid- 
ered in [21]. 
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FIG. 1. Analysis of the ground state at the critical point 
using the lowest Landau level approximation for A*' = 
12, g = 0.5 and A = 0.03. (Top) The coefficients 
Pn = \{N — 2n, 2n|\['Gs)|^ for the two-mode approximation to 
the many-body ground state at Oc ~ 0.776. The fidelity of the 
two-mode approximation is found to be |{^tm|*I'gs)|^ ~ 0.83. 
(Bottom) The quantum Fisher information as a function of 
rotation rate Q. 



When A ^ the anisotropic term connects subspaces 
of total angular momentum separated by two units of 
h. If the anisotropy is small, a finite maximum angular 
momentum Lmax can be chosen to ensure convergence 
for the energies and eigenstates of the the Hamiltonian 
restricted to the subspace of states with L < Lmax as 
considered in [19]. Throughout the rest of this paper, we 
use uj± as units of frequency and y/h/{Muj±) as units of 
length. 



III. RESULTS 

The case considered in [19] is that of a relatively slowly 
rotating gas at the threshold of the nucleation of the 
first vortex, where the strength of the interaction and the 
anisotropic perturbation are considered sufficiently small 
compared to the separation of the Landau levels in order 
to assume the LLL regime. In particular, the anisotropy 
strength is A = 0.03 {A <C 1) and the interaction one is 
assumed to be 5 iV = 6 for all the cases considered therein 
with different number of particles. This last expression is 
consistent with the standard mean-field criterion for the 
validity of the LLL approximation, /i <^ 2fiu}±. When 
Ng — 6 the typical rotational frequency for the nucle- 
ation of the first vortex is fJ > 0.77 and it can be shown 
that in unitless form [14], the standard criterion for the 
LLL validity is equivalent to Ng <C 27r/(l — fi), hence the 
Ng = 6 scaling is consistent with the LLL approximation. 
More importantly, it is also consistent with the power law 
scaling g^^y, = (6.92±0.08)iV-io46±o.oo5 ^j. g^^^jj ^^^_ 

ber of particles (TV < 20) describing the crossover from 
weakly to strongly interacting regimes for the first vortex 
nucleation as found by Feder [21]. 

As the rotation frequency is adiabatically swept, the 
system changes from being at rest to containing one vor- 
tex. In doing so, the ground state of the system passes 
a critical frequency f2c where the system undergoes tur- 
bulent symmetry breaking indicating a quantum phase 
transition. 

Moreover, for even numbers of particles, the approx- 
imation also predicts that at the critical frequency Jlci 
the ground state is a strongly correlated entangled state 
well-described by a two-mode approximation 
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where |A^ — 2n)|2n) is the state with {N — 2n) (respec- 
tively (2ri)) bosons in the the most (second most) popu- 
lated single-particle state ■i/'i (^'2)- These two states are 
eigenstates of the single-particle density matrix (SPDM) 

Pik ~ (^GslaJ^ailvE'Gs) and their populations are equal at 
the critical frequency, which together account for most 
of the population of the SPDM. Below the critical fre- 
quency, the most populated single-particle state ipi is ap- 
proximately a coherent superposition of two off-centred 
vortices with even parity corresponding to the states with 
m = units of angular momentum and to = 1 units. 
The second most populated state is a well-centred single- 
vortex state with odd parity corresponding to the state 
with m — 1. Above the critical frequency, ■01 and ip2 
abruptly swap places. 

The relevant results using the LLL approximation (as 
in [19]) for the current discussion are shown in Fig.l for 
A^ = 12 particles. The top plot shows the coefficients 
P„ = |C„|^ for the two-mode approximation of the many- 
body ground state in Eq.(3). This is promising for quan- 
tum metrology since the state has a "bat-like" structure 



which is known to be robust to particle loss in interfer- 
ometric schemes [18]. The bottom plot is the quantum 
Fisher information as a function of rotation rate f2. This 
is a measure of how well the quantum state can perform 
in metrology schemes. Suppose that the state given by 
Eq.(3) picks up an undetermined linear phase as a result 
of a unitary operation C/((^) being applied to the first 
mode, i.e. 
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The best precision, A0, with which (j) can be measured 
using this state in a single shot experiment, independent 
of the details of the measurement scheme, is given by the 
quantum Cramer- Rao bound [25], 
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where Fq is the quantum Fisher information. For a pure 
state 1^(0)) that depends on a single parameter (j), Fq is 
given by 

Fq = 4 [(Vl/'(0)|*'(0)) - |(VI/'(0)|*(0))|2] , (6) 

where 1*'(0)) = 9|*(0))/90. The results for the Fisher 
information in Fig.l show a relatively broad curve in the 
rotation frequency, which is a desirable property because 
it implies that experimentalists would have a sizeable 
margin of error when trying to hit the critical frequency 
to prepare the initial entangled state. 

Surprisingly, the form of the ground state and quantum 
Fisher information around criticality change dramatically 
when another Landau level is taken into account in the 
calculations. The results for the same set of parame- 
ters but using two Landau levels are shown in Fig. 2. Al- 
though the critical frequency has been shifted by roughly 
5%, the ground state at the critical frequency has now 
become more like a NOON state, (1^)10) + |0)|A^))/^/2, 
which can attain Heisenberg-limited precision that scales 
as A(f> ^ l/N but, unlike a "bat-like" state, is extremely 
fragile to particle loss. The fidelity of the ground state in 
the two-mode approximation with the full ground state 
is also reduced due to the increasing participation of a 
third single-particle mode which was not present in the 
LLL results. Furthermore, the quantum Fisher informa- 
tion shows an equally striking difference compared to the 
LLL approach. The Fisher information curve has become 
very sharp: its width has been reduced by about two or- 
ders of magnitude compared with Fig. 1. This sharp res- 
onance means that it is likely to be experimentally very 
challenging to produce a NOON state because a high con- 
trol of the rotation frequency is required in order to hit 
the critical frequency exactly. Here, it is important to 
mention that the quantum Fisher information has been 
calculated using the full ground state as obtained from 
the exact diagonalization procedure and after a single- 
particle change of basis to the natural orbitals one (for 
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FIG. 2. Analysis of the ground state at the critical point using 
a two Landau level approximation for A^ = 12, g = 0.5 and 
A = 0.03. (Top) The coefficients P„ = \{N - 2n,2n\^Gs)f 
for the two-mode approximation to the many-body ground 
state at Sic ~ 0.823. The fidelity of the two-mode approxi- 
mation is found to be |(*tm|*gs)P = 0.70. (Bottom) The 
quantum Fisher information as a function of rotation rate fl. 
Notice that the scale is roughly 10 times smaller than that of 
Fig. 1. 



which the SPDM is diagonal) , rather than using the two- 
mode approximation for the ground state. 

As it stands, these results do not look particularly 
promising for metrology. Fortunately, we can get all the 
benefits that the results in [19] suggest by looking more 
widely in parameter space. This is done in the next sec- 
tion where we assess the validity of the LLL approxima- 
tion in the current context. 



IV. VALIDITY OF THE LLL 

We have found in the previous section that, although 
the scaling A^^max = 6 provides a value for the maximum 
interaction strength gmax that is consistent with the use 
of the LLL approximation to accurately predict the crit- 
ical rotation frequency f2j, for a small number of atoms 
(A^ < 20), it is not a good guidehne to correctly pre- 
dict the quantum Fisher information or the exact form 
of the multi-body ground state at the critical frequency. 
This scaling relation is the result of thorough numeri- 
cal evidence for axisymmetric traps, where the criteria 



used to establish the validity of the LLL approximation 
is that the maximum interaction strength ^max predicts 
the rotation frequency fli for the nucleation of the first 
vortex within a relative error of 10% between the results 
of the LLL aproximation and those obtained with two or 
three Landau levels [21]. At this nucleation frequency ili, 
the ground state and first excited state become degener- 
ate and their energies show a crossing point. Therefore, 
in order to bring the condensate from the non-rotating 
ground state to the ground state containing one vortex, it 
is neccesary to break the rotational symmetry, thus cre- 
ating an avoided crossing that allows the passage of the 
system to the rotating regime with total angular momen- 
tum L — N and one vortex. The rotational symmetry 
can be broken using a weak (A <C 1) anisotropic stirring 
potential that helps to both impart angular momentum 
to the system and create the avoiding crossing. This ap- 
proach has been achieved experimentally for a few num- 
ber of atoms in optical lattices, where the atoms in each 
lattice are set in rotations using a local anisotropic stir- 
ring potential [24]. Naturally, it is reasonable to expect 
that the addition of this very small anisotropic perturba- 
tion (A <C 1) still allows for the description of the system 
in terms of the axisymmetric solutions. However, surpris- 
ingly enough, no matter how small A is (we checked this 
for values as small as oi A = 1 x 10~^, which is smaller 
than the typical static anisotropy defect in real traps), 
the striking differences in the quantum Fisher informa- 
tion and the exact form of the ground state found in the 
previous section still hold. Consequently, it is important 
to investigate the validity of the LLL approximation in 
the current context of the inclusion of small anisotropic 
perturbation. Some intuition as to the reason why the 
LLL approximation has failed to capture the physics pre- 
dicted by the addition of more Landau levels can be ob- 
tained by taking a closer look at the energy spectrum 
of the many-body Hamiltonian H. In the absence of 
the anisotropic term (A = 0), the interacting many-body 
Hamiltonian is exactly diagonalizable in blocks of definite 
angular momentum since the interaction term commutes 
with the total angular momentum (L^ in the present case 
of a 2D gas) and therefore does not mix subspaces of dif- 
ferent angular momentum, consequently, the eigenvectors 
of the interacting Hamiltonian are linear combinations 
of angular momentum eigenstates with the same eigen- 
value of Lz- In the LLL approximation, all the lowest 
eigenstates of H with L = 0, 2, . . . A^ become degenerate 
at the first energy crossing regardless of the interaction 
strength g and the ground state and first excited state 
consist of the states with L ~ and L = 2 respectively 
for rotation frequencies below and near the frequency of 
the crossing and the states L = N and L — N — 2 for 
frequencies above and near the frequency of the cross- 
ing. The effect of adding a small anisotropic term to H 
is that of lifting the degeneracy and nucleating an enta- 
gled "bat-like" state at the critical frequency fie which is 
mainly a superposition of the ground state L — N and 
excited state L = A^ — 2 of the unperturbed Hamilto- 



nian. These two states are directly connected by the per- 
turbation and they are energetically close together, thus 
making the superposition favorable. In contrast, when 
another Landau level is included, the isotropic eigen- 
states with L = 0,2,... A'' no longer cross at the same 
frequency, however, they remain quasi-degenerated for 
small values of g ~ (6/./V) x 0.4 and thus the addition 
of the anisotropic term has the same effect as in the 
LLL case, i.e. the nucleation of a "bat-like" state. In 
fact, the entangled ground state at f2c is again mainly 
a superposition of the states L — N and L = N — 2. 
On the other hand, when the interaction strength be- 
comes of the order oi g ^ (6/iV), the eigenstate L = in 
the isotropic case has a greater energy separation from 
the other L = 2 . . . N states which remain fairly quasi- 
degenerate at a certain frequency which is no longer the 
first energy crossing for the ground state, instead, the 
first crossing occurs above this frequency and consists 
only of the states with L — and L = N. In this case, 
the addition of the anisotropic term results in the nucle- 
ation of a "cat-like" state which is now mainly a super- 
position of the states L — and L = N. Although these 
states are not directly connected by A^ their proximity in 
energy and the fact that the states to which they are di- 
rectly connected through the anisotropy are farther away 
in energy, favors such superposition. 

In order to further quantify the validity of the LLL ap- 
proximation, we have considered the fractional change in 
two different quantities as given by the calculations per- 
formed with the LLL approximation and those obtained 
with more Landau levels included for A^ = 4, 6, 8, 10 and 
12 particles. The results are qualitatively similar for all 
the studied number of particles and we show the ones 
for TV = 12 for definiteness. the ones for TV = 12 for 
definiteness. First, we consider the predicted critical fre- 
quency flc- Our results shown in Fig. 3 are consistent 
with the scaling iV^max = 6 for axisymmetric traps. This 
critical rotational frequency can be predicted within a 
relative error of less than 10% when using the LLL ap- 
proximation and the two level approximation for a broad 
range of values of the interaction strength. Moreover, an 
extremely accurate prediction of the critical frequency 
is given by a two level approximation compared to the 
one obtained with three Landau levels. This picture is 
completely similar to the one that happens in axisym- 
metric traps when the fractional change in fii is con- 
sidered. A totally different picture is observed for the 
fractional change of the quantum Fisher information in 
Fig. 4. Here, we show the ratio (solid diamonds) be- 
tween AFq = \FQ{g)2hL - -Fq(5)lll|, the difference in 
Fisher information obtained with the LLL approxima- 
tion and that obtained with two levels at the critical 
frequency ilc obtained with a LLL calculation, and 
the Fisher information Fq obtained using a two Lan- 
dau levels approximation at this same critical frequency. 
Also, we show the same ratio (empty diamonds) for the 
case involving two Landau levels and three Landau lev- 
els. These results show vast disagreement in the Fisher 
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FIG. 3. The fractional change in the position of the critical 
rotation frequency Qc as a function of the interaction 
strength, calculated using two different approximations 
differing by one Landau level. In each case, the fractional 
change is calculated at the critical frequency obtained using 
the lower number of Landau levels. One level approximation 
(LLL) compared against the two level approximation (solid 

box) where AQc/^c = \i^c{g)2LL - i^c{g)LLL\/iic{g)LLL- 

Two level approximation compared against the 
three level approximation (empty box) where 

AD,c/^c = \flc{g)3-LL - f^c(ff)2LL|/f^c(ff)2LL. We Only 

show results for g > 0.2 since below this value, the partic- 
ipation of a rapidly increasing third single-particle mode 
deteriorates the fidelity of the ground state in the two-mode 
approximation and shifts the critical frequency above the 
trapping frequency which is physically unattainable. 
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FIG. 4. The fractional change in the quantum Fisher 
information Fq as a function of the interaction strength, 
calculated using two different approximations differing by 
one Landau level. In each case, the fractional change is 
calculated at the critical frequency obtained using the lower 
number of Landau levels. One level approximation (LLL) 
compared against the two level approximation (solid dia- 
mond) where AFq/Fq = \FQ{g)2LL ~ FQ{g)LhL\/FQ{g)2hL■ 
Two level approximation compared against the 
three level approximation (empty box) where 

AFq/Fq = |i^Q(sr)3LL - FQ{g)2LL\/FQ{gyiLL. 
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FIG. 6. The quantum Fisher information of the ground state 
for various basis truncations as a function of the interaction 
strength at the critical frequency Jlc calculated with the re- 
spective truncated basis. Three different basis truncation are 
considered: one Landau level (blue solid box), two Landau 
levels (red solid diamond), and three Landau levels (green 
solid circle). 



FIG. 5. Analysis of the ground state at the critical point using 
a two Landau level approximation for an interaction strength 
oi g = 0.2 and N — 12, and A = 0.03, which has been tuned 
for better suitability of the ground state for quantum metrol- 
ogy. (Top) The coefficients P„ ^ \{N - 2n, 2n|*Gs)P for the 
exact ground state at il^ ~ 0.938 calculated with the two level 
approximation. The fidelity of the two-mode approximation 
is found to be J(*tm|*gs)P = 0.80. (Bottom) The quantum 
Fisher information as a function of rotation rate fl for the 
LLL approximation (blue empty box) and two Landau levels 
(red empty box). 



information calculated with the LLL approximation with 
respect to that one obtained with higher Landau levels, 
which in principle are "closer" to the true value of the 
quantum Fisher information at that rotation frequency. 
This disagreement even holds for values of the interaction 
strength g which are typically regarded as pertaining to 
the weakly interacting regime for which the LLL approx- 
imation is considered valid. 

In particular, it is worth noting that for the value 
g — 6/N, i.e. g — 0.5 for the results in Figs. 3 and 4, the 
LLL approximation predicts a critical rotation frequency 
with a relative error of less than 6% with respect to the 
rotation frequency calculated with two Landau levels, 
however, the LLL prediction of the quantum Fisher in- 
formation has a relative error of more than 8000%. This 
is directly related to the fact that for these values of the 
interaction strength, the quantum Fisher information ob- 
tained with two Landau levels is a very sharp resonance 
as a function of g, whereas that obtained with LLL ap- 
proximation remains approximately broad across these 
values of the interaction strength. 



Notably, as the interaction strength is decreased much 
below the value oi g = Q/N, we recover the lowest Lan- 
dau level regime again. In particular we observe that 
for an interaction strength of about g ~ 0.2 in the two 
Landau level approximation, we get very similar results 
to those originally considered using the LLL approxima- 
tion in [19], i.e. a "bat-like" state with a broad quantum 
Fisher information curve. These results are shown in 
Fig. 5. This "bat-like" state is extremely interesting from 
the point of view of quantum metrology due to three main 
advantageous characteristics among other ones. First, 
the phase precision A0 scales as \/2/ ^JN[N + 2) and 
thus, for large number of particles it is comparable to the 
precision scaling of a "cat" state y/2/N which is known 
as Heisenberg-limited precision. This sub-shot-noise lim- 
ited precision for the "bat-like" state is already present 
even for modest number of particles N ~ 10. Second, it 
has been shown [18] that in the more realistic scenario 
where we have particle losses, the "bat" state retains a 
much higher sensitivity than the "cat" state, making it 
a robust choice in interferometric schemes. Finally, from 
the experimental point of view, the width of the Fisher 
information curve for the bat-like state is broader in the 
rotation rate J7. This is very desirable since it allows for 
a greater margin of error when trying to nucleate this 
entangled ground state. 

These findings are exciting for metrology, not only be- 
cause we can still use this system to create a robust 
"bat-like" state but we could also tune from "NOON-like" 
states to "bat-like" ones by only changing the interaction 
strength, which could be achieved in experiments using 
Feshbach resonances [23]. Tuning the system in order to 
get a "NOON-like" state can be challenging due to the fact 



that the Fisher information is a very sharp resonance for 
this case and our predictions for the critical frequency flc 
using two levels cannot exactly pinpoint the true critical 
frequency. However, we can always estimate a small fre- 
quency window where we expect to find the true critical 
frequency by calculating 51c using three and four levels 
and extrapolating our results to predict the many-level 
result. The quantum Fisher information at the calculated 
critical frequency for each approximation is shown in Fig. 
6. These results show that the two level approximation is 
expected to give very accurate predictions for the Fisher 
information of the full many-body ground state at the 
true critical frequency. 

According to our simulations for small number of par- 
ticles (N=4. . . 12) using two Landau levels, the width of 
the Fisher information for a fixed number of particles as 
calculated using the left width of the curve at half the 
maximum decreases with increasing interaction strength 
in an exponential manner, and follows approximately the 
same curve for each number of particles considered in our 
simulations, going typically from broad value of 5 x 10^"^ 
for g = (6/iV) x 0.4 ("Bat-like") to value of the order 
of 5 X 10-5 for g = {6/N) ("Cat-like"). Similarly, for a 
fixed value of g, the width also decreases exponentially 
as the number of particles is increased but the gradient 
is dependent on the exact value for g. 

The quantum Fisher information at the calculated crit- 
ical frequency for each approximation is shown in Fig. 6. 
These results show that the two level approximation is 
expected to give very accurate predictions for the Fisher 
information of the full many-body ground state at the 
true critical frequency. As Fig. 6 already suggests, we 
have found that the form of the entangled state and the 
shape of the quantum Fisher information do not change 
substantially when going from the two-level approxima- 
tion to a three-level one for all values of the interaction 
strength considered in our calculations. In other words. 



when the shape of the entagled state given by a two- 
level calculation is "bat" -like (broad Fisher information 
curve), it remains "bat" -like using three levels and the 
converse is also true for a "cat" -like shape (sharp Fisher 
information curve). 



V. CONCLUSIONS 

Rotating matter waves in anisoptropic potentials could 
prove to be an exciting system for generating useful en- 
tangled states for metrology schemes and studying quan- 
tum phase transitions. However, caution needs to be ex- 
ercised when choosing experimental parameters. While 
theoretical studies show that the LLL approximation pre- 
dicts the critical frequency well, a more in-depth analysis 
shows that it predicts other features very poorly. Indeed 
the bat-like state predicted at the critical frequency by 
LLL is seen to be more NOON-like and has a sharp reso- 
nance in its Fisher information meaning that it is unlikely 
to be practical for metrology. However we have been 
able to identify another parameter regime that gener- 
ates an entangled state with the features we would want. 
This system has a rich structure and may prove to be a 
promising tool for engineering a range of useful entan- 
gled states. The next step is to consider the dynamics 
of the system and a detailed scheme for demonstrating a 
proof-of-principle quantum-limited gyroscope. 
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